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Abstract 

In this paper we study a two-photon time-dependent Jaynes-Cummings 
model interacting with a Kerr-like medium. We assumed that the electro- 
magnetic field is in different states such as coherent, squeezed vacuum and 
pair coherent, and that the atom is initially in the excited state. We studied 
the temporal evolution of the population of the excited level, and the second 
order coherence function. The results obtained show that this system has 
some similarities with the two-mode Stark system. We analize two photon 
entanglement for different initial conditions. 
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I. INTRODUCTION 

One of the most fundamental models in quantum optics is the interaction of a single two- 
level atom with a single quantized mode of radiation, described by the Jaynes-Cummings 



Hamiltonian Despite being simple enough to be analytically soluble in the rotating- 

wave approximation, this model is able to describe the quantum-mechanical aspects of the 
interaction between light and matter. It has led to nontrivial predictions, such as the 
existence of "collapses" and "revivals" in the atomic excitation ||, and has also allowed 
a deeper understanding of the dynamical entangling and disentangling of the atom-field 
system in the course of time Q. The interest in the Jaynes-Cummings model (JCM) is 
mainly due to the fact that some of its predictions are nowadays accessible to experimental 
verification ||. A JCM interaction can be experimentally realized in cavity-QED setups 
and, as an effective interaction, in laser-cooled trapped ions. 

During the last decade many theoretical and experimental efforts have been done in 
order to study two-photon processes involving atoms inside a cavity, stimulated by the 
experimental realization of a two-photon cascade micromaser || Two-photon processes are 
also an efficient way of generating squeezed states of the electromagnetic field. It has also 
been established that two-photon degenerate atoms inside a Kerr-like medium can generate 
squeezing amplification 0. The recent discovery of new materials with very high Kerr 
coupling H opened the possibility for the implementation of new experiments that generates 
entangled states [[| that can be used for a perfect Bell-state discrimination. 

A system composed by a three level atom in S, A, and V configurations interacting with 
two modes of the electromagnetic field was proposed and studied more than ten years ago by 
Yoo and Eberly [|Kj. Following these guidelines, Gou [|ll|] investigated the H configuration 



when the intermediate level can be adiabatically eliminated [0 . This approximation turns 
the original bilinear photon-level interaction into a three-linear one (usually called "non- 
linear non-degenerate two-photon" interaction). This model has been broadly used in order 
to study the time evolution of the atomic and photon operators, the second order coherence 
function, the one and two-mode squeezing, the atomic-dipole squeezing and the emission 



spectra |ll|,p!3HT8| . Usually, the two- level system has been considered initially in the excited 
state and the two-photons have been chosen initially in two independent coherent states 
TT| , [T5|JT6| , p!8[| , two mode squeezed states fill , pair-coherent states [jl3| , correlated SU(1,1) 



coherent states fll4| , etc. The A configuration, also called the Raman coupled model, when 
the intermediate level can be adiabatically eliminated was studied by Abdalla, Ahmed and 
Obada [B3], and independently by Gerry and Eberly p0| . There have been investigations 
of the atomic inversion, the appearance of antibunched light, the violations of the Cauchy- 



Schwartz inequality, population trapping, and squeezing jT9]-|23[]. Some similarities in the 
Rabi frequencies of both configurations for some special conditions of the parameters were 
reported f[4| , and the connection between these configurations was deeply studied in Ref. 



|24fl . It was reported that for the H [25| , ^6 and A [27| , |28| models, when doing the adiabatic 



approximation, the appearance of Stark shifts must be taken into account. In Ref. |24[ the 
Stark shifts were neglected in order to separate their contribution to the nontrivial dynamics 
studied. The intensity-dependent Stark shifts caused by off-resonant levels J29| were modeled 
by Moya-Cessa, Buzek and Knight 0] using a JCM with an intensity-dependent shift of 
the two-level transition. Recently, some similarities and differences between the models 
used in order to describe a cavity filled with a Kerr-like medium, modeled by an anharmonic 
oscillator, (first analyzed in detail by Jex and Buzek pi| and discussed in many articles 
thereafter ) and Stark effects have been studied 

In this paper we study a two-photon time-dependent Jaynes-Cummings model interacting 
with a nonlinearity that can be, for instance, a Kerr-like medium. The problem is solved 
using a technique based on obtaining those observables dynamically related p4j, and then 
solving their temporal evolution. We assumed that the electromagnetic field is in different 
states such as coherent, squeezed vacuum and pair coherent, and that the atom is initially 
in the excited state. We studied the temporal evolution of the population of the excited 
level, and the second order coherence function. The results obtained show that this system 
has some similarities with the two-mode Stark system. 
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II. THE SYSTEM 



We study an effective two-level atom [24], whose levels \g > and |e >, with energies E\ and 
E 2 , interact with two modes of electromagnetic radiation of frequencies u)\ and uo 2 inside 
a non-linear Kerr-like medium. In two-photon processes there are more than two levels 
involved, but it is possible to neglect them if we assume that ui\ -\- uj 2 ~ E\ — E 2 (h = 1) and 
we consider that the transition frequencies between |e >, \g > and the intermediate levels are 
different from the frequencies of the field. So, the JCM in the rotating wave approximation, 
(H configuration), reads P4] : 



H = ^Eib% + ^tUidjai + T(t) (7010361 b\ + ^*b 2 b\a 2 d\) , (1) 

i=l i=l 

where a^, a\, b\ and 6«, i = 1,2, are creation and annihilation bosonic and fermionic operators, 
respectively, and 7 is the coupling constant between the atomic levels and the fields. 

The Kerr medium can be modeled by an anharmonic oscillator coupled to the field 
3T H^p9jl^[37f • Using the adiabatic approximation {u) me di um <^ to field) the non-linear 



medium can be represented by a non-linear power of the field. Thus, the total Hamilto- 
nian reads: 

2 2 



H = y^ j E i b\b i + V^ajai + T(t)(7aia 2 6i& 2 + Yb^lalal) 

i=l i=l 

+Xi(a t i) 2 (ai) 2 + X2(al) 2 (a 2 ) 2 + 2^x^ 2 a\a\a x a 2 (2) 



It is important to notice that three non-linear terms appear due to the presence of the non- 
linear Kerr media. Two of them are similar to the ones appearing in the case of one mode 
while the third one is a bilinear connection between the modes. This bilinear interaction is 
the one proposed in Refs. as a way of generating cat's states. 

The dynamical evolution of the mean values of the operators can be obtained by using 
the Ehrenfest equation, 



d(Pj)t 
dt 



5>^>, (3) 



where gij are the structure constants of a semi-Lie Algebra closed under commutation with 



the Hamiltonian [23], i.e. 

H{t)A] =ihY,9ji(t)d j . (4) 

i 

The operators defined via the previous equation are called relevant operators (RO) |]38| . The 
relevant operators for the two-modes Jaynes-Cummings Hamiltonian read 

iVT' m = (a t 1 )"(4) m 6 T 1 6 1 (a 2 ) m (a 1 ) n (5a) 

^ m = (air^rbik^r^r , (5b) 

A?' m = (a t 1 ) n (a t 2 ) m a t 1 a 1 (a 2 ) m (a 1 ) n , (5c) 

A? m = (alnalrata^riaxr , (5d) 

j„, m _ (at)»(4)"'( 7 a 1 a 2 Si^ + 7*44^1^) (a 2 ) m (ai) n , (5e) 

F n ' m = (a t 1 ) n (4) m z( 7 aia 2 S 1 SS - 7*4a 2 6 t 1 6 2 )(a 2 ) m (a 1 ) n , (5f) 

JV£™ = (a t 1 )"(4) m 6^ 2 5t^ i(a2)m(ai) n (5g) 

Notice that J n ' m and _F n ' m are the operators that have the information about the correlation 
between the modes and the atom levels (i.e. the entanglement). 

Using the Ehrenfest theorem, the evolution equations of the relevant operators are 

d (ivr m \ 

1 -T(t)F n ' m + nT(t)F n ~ 1 > m + mT(t)F n ' m - 1 

+nmT{t)F n - 1 ' m - 1 , (6a) 

-T{t)F n ' m , (6b) 
(n + l)T{t)F n ' m + mT{t)F n ' m - 1 

+mT(t)F n+1 ' m - 1 + nmT(t)F n > m - 1 , (6c) 
(m + l)T(t)F n ' m + nT(t)F n - 1 ' m 

+nT{t)F n - 1 ' m+1 + nmT{t)F n - l ' m , (6d) 

, (6e) 



dt 



d 



(iV 2 "< m ) 
dt 



dt 



d(A^ m 



dt 



d(N^ 
dt 



d ( I n '' 



dt 



d ( F n '' 



dt 



[a -2 {n X i +m X 2 + (m + n + l)Vxi*i)] i> ' m 
-2 (Xi + VXVS) F n+1 ' m - 2 (X2 + VXV& F n ' m+1 , 

- [a - 2 (n X i +m X 2 + (m + n + l) ^xl)} I n ' m 

+2 (xi + Vxljft) ^ +1 ' m + 2 (X2 + VxDfc) ^"' m+1 
+2| 7 | 2 T(t)[(n + l)(m + l)iV 2 n ' m - (n + l)(m + l)iV£f 
+(n + l)iV 2 n ' m+1 - (n + l)iV£f +1 + (m + l)iV 2 " +1 ' m 



(6f) 



(m + 1)KJ 



n+l,m 



jy-n+l,m+l _|_ jy-n+l,m+lj 



(6g) 



where a = E 2 — E\ — u\ — uo 2 . 

So, using Eqs. @ it is possible to evaluate the temporal evolution of the RO. These 
operators close a Lie Algebra. Thus, any extension of the two-photon JCM built up adding 
terms proportional to RO will have similar evolution equations. This fact will be used in 
the following section in order to study the two-photon JCM with Stark shifts. 

Equations (|]) can be numerically solved for any temporal dependence of the atomic-field 
interaction, T(t). For the time-independent case, Eqs. (§) can be analytically solved using 
the series expansion in terms of commutators with the Hamiltonian p^j3^ j3 

1 / t x " 



O 



6 +E 



o * — ' n\ V ih 

n>0 



0,H 



H 



(7) 



From Eqs. (0) we can notice that the evolution equations for all the operators different from 



pn,m depend on pn,m^ g Q ^ double commutator 
double commutator reads 



H, 



H,F r 



plays a central role. This 



H,F r 



_ f>2 pn,m i (o2 _ o2 \ pn+l,m , 

r J n,m 1 \rn+X,m t^n,mj 1 

(Plm+i ~ Pl,m) F n ' m+1 + [8e l£2 + 4| 7 | 2 ] F n+1 > m+1 + 

where the generalized Rabi frequency /3 2 m = [a — 2 {n\i + mX2 + (n + m + 1)^X1X2)] 2 + 
4|7| 2 (n + l)(m + l) and e« = x% + \fXXXi- ^ is important to notice that the term propor- 
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tional to ^/XiX2 appearing in 0^ m and is due to the bilinear connection between modes. 
Equation @ generates paths, or different ways of dynamically connecting the RO, repre- 
senting quantum correlations pi,51| S5|,P5| . From Eq. (|8|) it can be seen that those terms 
proportional to q are due to the presence of the Kerr medium. Notice that taking into 
account the bilinear connection between modes doubles the value of e« in the xi — X2 case 
usually studied. 

Making use of Eq. (H) we can obtain the temporal evolution of the population of 
the excited level ((A^' )^ and the second order intermodes coherence function (defined 
as gj 2 (t) = ,^ 1 , ai J 2 t a2 \ — 1), for different states of the field, (See appendix). 

This is depicted in Figs. [l]-[3]. In all cases, we assume the atom initially in the excited state 
((A r { ) ' )(0) = 0, (iV2' )(0) = 1). The phenomenological third-oder nonlinear susceptibility 
for the Kerr media X = Xi = X2 takes the values 0, 0.5 and 1, which, in principle, could 
be obtained in the ultra-high susceptibility Kerr cells recently studied ||. We also assume 
that there are initially 10 photons in each mode (|o!j| 2 = 10). We analyze the influence of 
the non-linear medium in the temporal evolution of the excited level and we find inhibited 
decay as it can be seen from Figs. [I|-b,^|-b,Q-b. In the Pair Coherent state case (Fig. |3|-b), 
we find that the revivals are of such a kind that the atom recovers its initial population 
in each revival, and this effect is increased as we make more important the coupling with 
the non-linear medium. Another interesting feature is that in the squeezed vacuum state 
(Fig. 0-b) we see that the revivals are regular and sharp, and become more periodic as we 
increase the strength of the coupling with the Kerr medium. We can see from Figs. |l|-a and 
0-a, that the two-photon JCM presents antibunching, and that this effect is attenuated by 
a stronger interaction with the non-linear medium. In the particular case of the squeezed 
vacuum, we notice that the field recovers its initial intermodes coherence when the excited 
level population reaches a maximum. 

Finally, we want to mention that A configuration can be studied by using the canonical 
transformation di = a\ and redefining the RO following Ref. 
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III. STARK EFFECT 



The Stark effect is related to changes in the atomic energy levels due to virtual transitions 
from levels out of resonance. The effective Hamiltonian taking into account the Stark effect 
reads 

2 2 



i=l i=l 

+a\ai ypib^i + 772^2^2 j + a\a 2 {vib[bi + rj 2 b\b 2 ^ . (9) 
Both the Kerr Hamiltonian and the Stark Hamiltonian can be written in terms of RO (for 



the one-mode case see Ref. [pq| . Using the same technique is is easy to see that if 



^1 + VX1X2 = u si + m 



^2 + VX1X2 = w S2 + 772 

Vl -r] 2 = 2ei = 2e 2 , (10) 

the evolution equations for the relevant operators will be the same. In this case, for the RO 



here studied, both systems will have the same temporal evolution This will not be the 
case if we study the temporal evolution of other R.O. like the electromagnetic field (a + a^). 



IV. CONCLUSION 

In the present paper we analyzed the non-degenerate two-photon JCM in the presence 
of a Kerr media. We have identified those relevant operators that are interrelated with the 
dynamics of the system. We find that the dynamical evolution of the relevant operators 
for this problem and for the model that takes into account Stark shifts are the same for 
these RO and for some values of the constants. We analyze the influence of the non-linear 
medium in the temporal evolution of the excited level and we find inhibited decay. In the 
Pair Coherent state case, we observe that the revivals are of such a kind that the modes 
of the field are disentangled, and this effect is increased as we make more important the 



S 



coupling with the non-linear medium. We also notice that in the squeezed vacuum state 
the revivals are regular and sharp, and that the modes of the field are entangled. We also 
see that this system presents antibunching, and that this effect is attenuated by a stronger 
interaction with the non-linear medium. Thus, this system can be used in order to obtain 
very different two photon entanglement situations, depending on the initial state of the field 
that is considered. This fact can be very useful in experimental realizations. 
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APPENDIX: EXACT EXPRESSIONS FOR THE TIME EVOLUTION OF THE 

RELEVANT OPERATORS 

In the time independent case, making use of equations (||[7|) we can obtain an exact 
expression for all the Relevant Operators. For instance (A r { ) ' )(t) reads: 

OC ~l CO CO 

<jv?°)(t) = (ivD(o) -E^EE «;,f(^' m )o 

j,k=0 ^ n=j m=k 
OO ^ i oo oo 

- E 7# E £ <n 2ei + ^ - [a - 2(n £l + me 2 

j,k=0 3tk n=j m=k 

+ V^)W n ' m )o + 2| 7 | 2 [(ivr 1 - m+1 > - <iv7 +1 ' m+1 >] 

+ 2\ 1 \\n + l){m+l)[(N^ m )-(N^)} 

+ 2| 7 | 2 (n + i)[(iVT m+1 > - (K;r +1 )] , (Al) 

where S Jt k = sin((3 j>k t), C j>h = cos((3 jtk t)-l, and ojf = (-l) n+m +^ k+1 /[(n-j)\j\(m-k)\k\}. 
We can evaluate this expression for different initial conditions for the population of the 
atomic levels, and the state of the field. 

If we assume that the field is initially in a coherent state, 
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where 



\a x a 2 >= } J Pn 1 Pn 2 \nin2 > 



8=0 



el°l 2 b| 2ni 



then the expression (\KTj ) takes the form: 



1 C, 



(iVD, = 1 + 2| 7 | 2 ^ -^^(Ar) J o(Ar)Se-^ 0)o e-^°)o 



j,k=0 ■ 



where A, 0,0 (i = 1, 2) is the number of photons present in each mode of the field. 
The same treatment can be applied to the case of having the field in the 



Vacuum state p9fl , which can be written in the following way, 



",<f> >= (cosh(r)) 1 '^2(e l<l 'tanh(r)) n \n,n > . 



In this case the expression (|A1|) reads, 



(iVf °> t = 1 + 



2|7l 



cosh 2 (r] 



^(k + lftanh 2 



k=0 



Ck,k 

yk,k 



Finally, we assume the field is in the Pair Coherent state j|0 



\£, q >=N q J2 



n=0 [n!(7i + g)!)]3 



-\n + q,n > 



where N q is a normalization constant. In this case we have: 



(JVf°> t = 1 + 2| 7 | 2 iV g + l)(k + q + i)_|L_C7 M 



fc=0 
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